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Abstract. Recently it has been shown that under mild conditions Voronoi dia- 
grams have a certain continuity property: small perturbations of the sites yield 
small perturbations in the shapes of the corresponding Voronoi cells. However, an 
important assumption needed for formulating this result is that the ambient normed 
space is uniformly convex. Unfortunately, simple counterexamples show that if the 
uniform convexity is removed, then instability can occur. Since Voronoi diagrams 
in normed spaces which arc not uniformly convex do appear in theory and practice, 
e.g., in the plane with the Manhattan distance, it is natural to ask whether the 
stability property can be generalized to them, perhaps under additional assumptions. 
This paper shows that this is indeed the case assuming the unit sphere of the space 
has a certain (non-exotic) structure and the sites satisfy a certain "general position" 
condition related to the unit sphere. 
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1. Introduction 

1.1. Background. Given a world (space) X, together with a distance function and a 
collection of sites (subsets) {Pk)keK in the world, the Voronoi cell associated with 
the site is the set of all the points in X whose distance to is not greater than their 
distance to the other sites Pj,j ^ k. Voronoi diagrams have many theoretical and 
practical applications in various fields [7, 8, 45, 22] and therefore have been investigated 
quite a lot. However, most of this investigation has been focused in finite dimensional 
Euclidean spaces with point sites, and in many cases only in and M^. Research 
works studying these diagrams in spaces which are not Euclidean do exist, e.g. [4, 
10, 11, 17, 18, 30, 34, 35, 36, 37, 40, 49, 52, 55], but they constitute a small fraction 
comparing to the research in the Euclidean case, and, in fact, the majority of them 
study the diagrams in algorithmic aspects. As a result, not much is known about 
Voronoi diagrams beyond Euclidean spaces. 

This paper studies Voronoi diagrams in a class of finite dimensional normed spaces 
which are not necessarily Euclidean. More specifically, it studies a certain continuity 
property of the cells, namely that small changes of the sites, e.g., of their position or 
shape, yield small changes in the shapes of the corresponding Voronoi cells. Recently 
[50] it was shown that the Voronoi cells do have such a stability property, and it 
actually holds in a quite general setting (e.g., infinitely many sites of a general form 
in a possibly infinite dimensional space). However, an important assumption needed 
for formulating this result is that the space satisfies a certain property called uniform 
convexity. Roughly speaking, it means that not only the unit sphere does not contain 
any line segment, but actually that given any line segment whose endpoints are on 
the unit sphere, if the length of this segment is bounded below by a certain positive 
number e, then the middle point of this segment must penetrate the unit ball by at 
least a positive number 5(e), uniformly for all such line segments. While, in particular, 
the familiar Lebesgue spaces Lp{fl) and sequence spaces ip, 1 < p < oo are uniformly 
convex [19, 21, 39] , there are useful and not less familiar spaces which are not uniformly 
convex. Among them, (m > 2) with the £i (Manhattan) and i^o (max) norm. 

Voronoi diagrams based on non-uniformly convex norms (mainly ii and ioo) do oc- 
cur in theory and practice. A simple example is for giving a rough estimate on the 
number of costumers of a given facility (e.g., a shopping center/post office) located, 
together with its competitors, in a fiat region in which transportation is restricted 
to lines parallel to the standard axes. This is the case of Manhattan. See Figures 
1-2 for an illustration. A few additional examples can be given, e.g., in computer 
graphics [27] in which the sites are certain points in an image and the (centroidal ii) 
Voronoi cells are tiles having good shapes which are the building blocks for a mosaic; 
for storage [38] (•^i,-^oo); in relation with VLSI design [47, 48] in which the sites are 
non-point components and the i^o Voronoi cells allow one to estimate critical areas 
in the circuit; in databases [46], where the sites are vectors in a vector space and the 
intersections of the Voronoi cells with respect to varying ip norm p — 1, 2, 3, . . . , oo 



THE GEOMETRIC STABILITY OF VORONOI DIAGRAMS 



3 




Figure 1. 10 shops/post offices in a 
flat city, modeled as points. Distance 
are measured using the ii distance. 
The Voronoi cells of the sites are dis- 
played. 




Figure 2. A more realistic descrip- 
tion: the sites have slightly moved 
(since their location is not known ex- 
actly) and their shapes are not an ideal 
point. But the Voronoi cells have al- 
most the same shapes as before. 



are used for constructing a certain index structure; and possibly in analyzing (approx- 
imate) nearest neighbor search algorithms in high dimensional spaces with the ii or 
the ioo norms [3, 5, 29]. 

In all of the above examples imprecision (noise, input or computational errors, etc.) 
is an integral part of the setting, but despite this, analysis of the geometric stability 
of the cells seems to be absent here and in most of the literature dealing with Voronoi 
diagrams (with the exceptions of the very brief and intuitive discussions in [2], [7, p. 
366], [31] in Euclidean settings; see [50]). As a matter of fact, issues related to the 
stability/robustness of geometric structures and algorithms due to imprecision are 
not so common in the literature, at least comparing to the majority of works in which 
ideal conditions (infinite precision of the input, infinite precision in the calculations, 
etc.) are assumed. There are works which do deal in one way or another with such 
issues, e.g., [1, 2, 6, 9, 12, 13, 14, 20, 24, 25, 26, 28, 33, 41, 42, 51, 53, 54, 56, 57], but 
the focus is mainly on issues related to combinatorial properties (in Euclidean spaces). 

To the best of our knowledge, even when geometric stability is discussed in settings 
closely related to Voronoi diagrams, e.g., in [16] (medial axis) and focuses on the 
Hausdorff distance, as here, the question of stability of the shapes of the Voronoi cells 
under small perturbations of the sites have not been addressed or even raised at all. 
In addition, the results in such papers do not imply our ones even for the Euclidean 
case (for instance, in [16] one cannot take the set fl considered there to be the union 
of the sites or the complement of this union, because fl should be a bounded open set 
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with a boundary which is a smooth enough manifold, while usually neither the union 
nor its complement are such a set). 

Taking into account the above, one may be interested in having a theoretical result 
ensuring the geometric stability of the Voronoi cells in normed spaces which are not 
uniformly convex. At first it seems that establishing such a result is impossible, since 
there are simple counterexamples even in the plane with the Coo norm showing that in- 
stability of the Voronoi cells may occur under arbitrary small perturbations of the sites 
[50] (see also Example 6.1). However, a careful inspection of these counterexamples 
shows that the sites form certain degenerate configurations. Therefore, it is natural 
to ask whether in the common case, i.e., whenever such degenerate configurations are 
not formed, a stability property still holds. 

1.2. Contribution of this paper: We consider Voronoi diagrams in a class of finite 
dimensional normed spaces having a certain non-exotic unit sphere and show that: if 
the finitely many compact sites satisfy a certain "general position" assumption related 
to the structure of the unit sphere, if there is a positive lower bound on the distance 
between them, and if the discussion is restricted to a compact and convex subset, then 
the Voronoi cells are stable with respect to small changes in their corresponding sites, 
where the changes are measured with respect to the Hausdorff distance. Roughly 
speaking, the property which characterizes the class of normed spaces we consider is 
that their unit sphere has the finite face decomposition: it has a (possibly empty) 
"rotund part" and finitely many (possibly 0) "flat parts" each of them is a nonde- 
generate closed and convex subset. In particular, (]R™,£i) and (^"^,£00) ai'c in this 
class. The general position assumption is that no two points from different sites form 
a line segment which is parallel to a line segment contained in the unit sphere. If 
the sites are points and they are generated using the uniform distribution, then this 
condition holds with high probability (probability 1, unless computational errors and 
the discrete nature of computational devices are taken into account). 

Comparing to the approaches discussed in the references mentioned a few para- 
graphs above, our approach is a physical-analytical one: we do not measure any 
combinatorial difference between the real and perturbed Voronoi cells as a function of 
some measure between the sites but rather we can ensure that if the Hausdorff distance 
between the real and the perturbed sites is small enough, then so is the Hausdorff dis- 
tance between the real and the perturbed Voronoi cells. Since however we do not 
give explicit bounds (in contrast with [50] where the bounds are even dimension free, 
namely no "curse of dimensionality" holds), this continuity result is more theoretical 
in nature, analogous with establishing the convergence of a sequence/algorithm to a 
limit, but without giving estimates on the convergence speed. It is somewhat resem- 
bles the continuity result of Groemer [23] (but significant differences exist: see [50] for 
a comparison) . 

1.3. The structure of the paper: In Section 2 basic definitions and notations are 
presented. In Section 3 the concept of finite face decomposition is introduced and 
discussed. The general position condition is discussed in Section 4. The stability 
theorem is presented in Section 5 and some aspects related to its proof are discussed. 
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In Section 6 a few relevant counterexamples are presented, showing the necessity of 
the assumptions imposed in the theorem. The proof of the stability theorem can be 
found in Section 7. We end the paper in Section 8 with several concluding remarks. 

2. Notation and basic definitions 

Throughout the paper, unless otherwise stated, X is a nonempty compact and 
convex subset oi X — W^, m > 1, with some norm | ■ |. The induced metric is 

d{x,y) = \x — y\. The unit sphere is Sj^ = {^gX:|^^| = 1}. We assume that X is 
not a singleton, for otherwise everything is trivial. We denote by [p, x] and [p, x) the 
closed and half open line segments connecting p and x, i.e., the sets {p + t{x — p) : 
t e [0, 1]} and {p + t{x —p) : t E [0, 1)} respectively. The segment \p,x] is called non- 
degenerate if its endpoints p and x are different. Any (real) line L can be represented 
as L = {p + t9 : t E M.} where 6 (the direction vector) is any non-vanishing vector 
which, after normalization, can be assumed to have norm 1, and p is some point on L. 
A hne L = {p + t9 : t e R} is parallel to a second hne L' = {p' + t'9' :t'eR} if it is 
a translation of it by some vector, i.e., L = q + L' — {q + v' : v' E L'} for some vector 
q. Elementary manipulations show that this is equivalent to either 6 = 6' or 6 = —6'. 
Two non-degenerate segments [p, a] and [p', a'] are said to be parallel if they are located 
on parallel lines, and this is equivalent to either (a — p)/\a — p\ = (a' — p')/\a' — p'\ or 
{a-p)/\a-p\^-ia'-p')/\a'-p'\. 

Definition 2.1. Given two nonempty subsets P,A<ZX, the dominance region dom(P, A) 
of P with respect to A is the set of all x & X whose distance to P is not greater than 
their distance to A, i.e., 

dom(P, A)^{xeX : d{x, P) < d{x, A)}. 

Here d{x, A) = mf{d{x, a) : a E A} and in general we denote d{Ai,A2) = inf{(i(ai, 02) : 
Qi e Ai, a2 e A2} for any nonempty subsets Ai,A2. 

Definition 2.2. Let K be a set of at least 2 elements (indices). Given a tuple {PkjkeK 
of nonempty subsets Pk C X , called the generators or the sites, the Voronoi diagram 
induced by this tuple is the tuple {Rk)keK such that for all k & K, 

Rk = dom(Pfc, IJP,) ^{xeX: d{x, Pk) < d{x, Pj) Vj e K, j k}. 

In other words, the Voronoi cell associated with the site P^ is the set of all x E X 
whose distance to P^ is not greater than their distance to the union of the other sites 

We finish this section with the definition of the Hausdorff distance. 

Definition 2.3. Given two nonempty sets Ai, A2 C X, the Hausdorff distance between 
them is defined by 

D(Ai, A2) = max{ sup ^(01,^2), sup d{a2,Ai)}. 
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Note that the Hausdorff distance D{Ai,A2) is definitely different from the usual 
distance 

d{Ai,A2) = inf{rf(ai, 02) : ai G Ai, a2 G A2}. 

As a matter of fact, D{Ai, A2) < e if and only if d{ai, A2) < e for any oi G Ai, and 
d{a2,Ai) < e for any 02 G A2. In addition, if D{Ai,A2) < e, then for any oi G Ai 
there exists 02 G A2 such that d(ai, 02) < e, and for any 62 G A2 there exists bi G 
such that d{b2, bi) < e. These properties can be used for giving an optical-geometrical 
explanation for the use of Hausdorff distance as a natural tool when discussing ap- 
proximation and stability in the context of sets: suppose that our resolution is at most 
r, i.e., we are not able to distinguish between two points whose distance is at most 
some given positive number r. If it is known that D{Ai,A2) < r, then wc cannot 
distinguish between the sets Ai and A2, at least not by inspections based only on 
distance measurements. As a result of the above discussion, the intuitive phrase "two 
sets have almost the same shape" can be formulated precisely: the Hausdorff distance 
between the sets is smaller than some given positive parameter (note that a set and a 
rigid transformation of it usually have different shapes) . 

3. Subsets with a finite face decomposition 

In this section we introduce and discuss the seemingly new notion of finite face de- 
composition, a notion which plays an important role in the stability theorem (Section 
5). 

Definition 3.1. The nonempty subset S C R"* is said to have a finite face decompo- 
sition if ^t^\Fi where t G {0, 1, 2,3,.. .} and: 

(a) Fi is a closed and convex subset of S for each i G {1, ■■■,£} and it is not a point; 

(b) for each i E {1, . . . ,i} the subset Fi has a certain maximality property: given any 
line segment [s, s'] C S, if [s, s"] C F^ for some s" G (s, s'], then the whole segment 
[s, s'] is contained in Fi; 

(c) The subset F^+i C S does not contain nondegenerate line segments and F^+iflFj = 
for any i G {!,...,£}. 

The subsets F^, i G {!,...,£} are called the flat parts and F^+i is called the rotund 
part. 

Example 3.2. The unit sphere of M"^ with the £^0 norm || (xi , . . . , ^c^) || — max-[ \xi | : 
i = 1, . . . ,m} is a cube. It can be decomposed into i = 2m rectangular faces. This 
shows that it has a finite face decomposition. Here F^+i = 0. By the same reasoning 
the unit sphere of R"* with the ii norm . . . = \xi\ has the finite 

face decomposition. The standard Euclidean sphere has a trivial finite face decom- 
position where i = and F^+i is the sphere itself (because it does not contain any 
nondegenerate line segment; so is the case in any strictly convex norm space such as 
the ip spaces, 1 < p < 00). In general, any m-dimensional compact polyhedron in 
R™ which is symmetric with respect to the origin and the origin is an interior point 
of it (with respect to the Euclidean norm) induces a new and equivalent norm on 
whose unit sphere is the boundary of this polyhedron. This unit sphere has the finite 
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Figure 3. Several 2D unit spheres Figure 4. The unit sphere described 
with the finite face decomposition: £i in Example 3.3. 

(rhombus), £2 (circle), £00 (square), 
something else (octagon). 

face decomposition with i=the number of faces, F^+i = 0. See Figure 3 for several 
illustrations. 

Example 3.3. Consider the unit sphere of the Euclidean cut by the planes X3 = a 
and X3 = —a. Let S = FiU F2U where a G (0, 1) is fixed and 

-^1 = {{xi, X2, X3) : xf + xl + xl < 1, X3 = a}, 
F2 = {{xi,X2, X3) : xl + xl + xl <1, X3 = -a}, 

F3 = {{xi, X2, X3) : xl + xl + xl = l, XsE (a, -a)}. 

See Figure 4 for an illustration. This is a finite face decomposition of S. Since S is 
symmetric with respect to the origin and the origin is an interior point of the body 
induces by S, it follows that 5* induces a new and equivalent norm on M^, and the 
corresponding unit sphere is 5*. This unit sphere has a finite face decomposition. 

Example 3.4. For each ri G N define a the pair (x„, i/n) G by x„ = 1/n and ?/„ = 
1— (1/n). Now connect by a line segment each of the pairs (a;„, yn) and (x„+i, Un+i)- Do 
the same with the pairs (— x„, —yn) and {—Xn+i, —yn+i)- Finally, connect (xocZ/oo) = 
(0,1) to {-xi,yi) = (-1,0) and also (xoo, -2/oo) = (0,-1) to {xi,yi) = (1,0). Let S 
be the union of all of these line segments. The set 5 is a polygon with infinitely many 
edges (all of them are contained in the first and the third quadrants) and its vertices 
are the pairs (±x„, ±?/„), G N U {00}. It induces a new unit sphere in the plane 
which does not have any finite set decomposition. Intuitively, this is because the sets 
Fi should be the line segments, but there are infinitely many line segments. 

Example 3.5. The "cylinder with covers" S = {{xi,X2,X3) : + |a;2p = 1, 0:3 G 
[—1,1]} does not have a finite face decomposition. Intuitively, this is because the 
cylinder contains infinitely many (vertical) line segments and each of them must be a 
face. 
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Figure 5. Illustration of Condition 
4.1: two sites in {R'^Joo)- 

4. The "general position" condition 

The goal of this short section is to discuss further the "general position" assumption 
the sites {Pk)k£K need to satisfy in order to allow the cells to be geometric stable with 
respect to small changes of the sites. The precise condition is stated as follows: 

Condition 4.1. For any j, k E K , j ^ k and for any Pj G Pj,Pk G Pk, Pj Pk, the 

nondegenerate line segment [pj,Pk] is not parallel to any nondegenerate line segment 
contained in the unit sphere of the space. 

An illustration of the condition in the case of with the £oo norm (and two sites) 
is given in Figure 5. To the best of our knowledge, Condition 4.1 has not been dealt 
with in the literature. However, in very few places, e.g., in [7, p. 390, Figure 37], [34, 
pp. 19-20], [37, p. 605, Fig. 1(b)], [45, p. 191, Figure 3.7.2], one can find hidden 
clues to its possible existence in the particular case of (M^, ii) and 2 point sites, either 
Pi = {(-1,-1|)} and P2 = {(1,1)} or Pi = {(-1,1|)} and P2 = {(1,-1)}. In 
this case the segment [^1,^2] generated by the (unique) points pi = (—1,-1) G Pi 
and P2 = (1, 1) G P2 is parallel to 2 segments contained in the unit sphere of the 
space (see Figure 3). Neither this fact nor any kind of discussion regarding a possible 
instability of the shapes of the Voronoi cells are mentioned in the references described 
above, but it is acknowledged there that the bisectors induced by this configuration 
are problematic and exotic (because they are "fat"). 

As a final remark regarding Condition 4.1, let us emphasize that the "general po- 
sition" property required by this condition is significantly different from the common 
"general position" condition frequently found in the computational geometry litera- 
ture dealing with Voronoi diagrams. In its simplest form, where the setting is the 
Euclidean plane, this condition says that no 3 distinct (point) sites are located on the 
same line and no 4 sites are located on the same circle. In contrast, in Condition 4.1 
there is no problem at all in both cases even if only point sites are considered. For 
instance, 3 point sites that are located on the same line will not induce instability as 
long as this line is not parallel to a line segment contained in the unit sphere of the 
space. 
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Figure 6. 10 point sites and their 
cells in a square in (M^,£oo)- Condition 
4.1 holds. 



Figure 7. The sites have slightly per- 
turbed. The shapes of the cells have 
slightly perturbed. 



5. The stability theorem and some aspects related to its proof 

In this section we formulate the stability theorem and discuss briefly issues related 
to its proof. For a simple illustration of the theorem, see Figures 6-7. 

Theorem 5.1. Let X be a compact and convex subset ofMJ^, m G N. Let \ ■ \ be a 

norm defined on the space and assume that the unit sphere of the space has a finite 
face decomposition. Let K be a finite set of indices and let {Pk)keK be a given finite 
tuple of nonempty compact subsets of X . For each k E K let Ak = [jj^k^J ^''^^ 
Rk = dom[Pk, Ak) be the Voronoi cell corresponding to Pk- Suppose that 

T] := mm{d{Pk, Pj) ■.j,ke K, ] ^ k} > 0. (1) 

Suppose also that the sites are located in a "general position" with respect to the unit 
sphere, namely, that Condition 4-i holds. Then for each e > there exists A > such 
that if {Pl^)keK is any tuple of nonempty compact subsets of X satisfying D{Pk,Pl.) < 
A for each k E K, then D{Rk, R'j.) < e for each k E K. 

The proof of Theorem 5.1 is quite long and technical, and it is given in Section 7. 
The proof itself is partly inspired by and partly based on arguments and assertions 
given in [50]. In particular, there is a heavy use of spherical arguments (unit vectors, 
unit sphere). However, there are considerable differences between both proofs and 
some of the involved ideas since a key step in [50], namely Lemma 5.5 in the extended 
abstract version and Lemma 8.9 in the current arXiv version, aimed at proving a 
certain geometric estimate, is heavily based on the assumption that the normed space 
is uniformly convex. Important tools used here in order to overcome the difficulty 
of lacking of this property are geometric arguments based on the compactness of 
the involved sets, arguments based on finiteness, a new geometric characterization of 
equality in the triangle inequality (instead of the forgotten strong triangle inequality 
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Figure 8. The figure of Example 6.1: 
20 point sites in a square in (]R^,£i). 
Only four cells are displayed. 



Figure 9. The non-center sites have 
slightly moved (in the first component) 
towards the center site. The verti- 
cal rays of the center sites have disap- 
peared. 



of Clarkson [19, Theorem 3] used in [50]), and a certain quantitative way of measuring 
how far a given vector is from being parallel to a non-degenerate segment contained in 
the unit sphere. The geometric characterization of equality in the triangle inequality 
mentioned above is as follows: 

Lemma 5.2. Let xi and X2 be two non-vanishing points in a normed space. Let 
Xi = Xi/\xi\, i = 1,2. Then \xi + X2I = + \x2\ if and only if the segment [xi, £2] is 
contained in the unit sphere. 

The proof of this lemma is based on another characterization of equality in the 
triangle inequality, observed recently by L. Maligranda [44, Corollary 1.5]: \xi +X2\ = 
\ + 1x2! if and only if Ixi + £21=2 (Maligranda's characterization is also a simple 
consequence of results mentioned in [32, 43]). 

6. Counterexamples 

In this section we discuss briefiy several counterexamples related to Theorem 5.1 
showing that the assumptions imposed in the theorem are necessary. 

Example 6.1. If the sites do not satisfy the "general position" assumption (Condition 
4.1), i.e., some of them form line segments parallel to segments contained in the unit 
sphere, then the Voronoi cells may not be stable as shown in Figures 8-9. Here 
the setting is with the ii norm. There are 20 point sites in the rectangle X = 
[—10,10] X [—10,10]. The sites form a symmetric structure composed of 4 groups 
of 5 sites in each group. The bottom left group consists of Pi = {(— 6, — 6)},P2 = 
{(-2,-6)},P3 = {(-6,-2)},P4 = {(-2,-2)} and the center site P5 = {(-4,-4)}; 
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the other sites are obtained by translating this group with the vectors (8,0), (0,8), 
(8,8). Only the cells of P5, Pio, P15, P20 are displayed. A small perturbation of the 
other sites so that they will be closer in the first component to the center site of 
the group by some arbitrary small /? > causes the vertical "rays" of these sites to 
disappear, but it preserves the horizontal rays (these rays are, actually, more or less 
very thin strips). 

Example 6.2. The positive lower bound expressed in (1) is necessary even in a 
square in the plane. Consider X = [—10,10]^, the ii distance, Pi^^ = {(0,/3)} and 
-P2,/3 = {(0, -/3)}, where /3 G [0, 1]. As long as /3 > 0, the cell dom(Pi,^, P2,/}) is the 
upper half of X. However, if /3 = 0, then dom(Pi^o, -^2,0) is X. This example is not 
specific to non-uniformly convex spaces and can be formulated in the Euclidean plane 
too, as actually done in [50]. 

Example 6.3. As for the compactness assumption, let X be the (non-compact) plane 
with the £00 norm and let Pi = {(0, lOn) : n G N}, P2 = {(l/n,5 + lOn) : n G N}. 
Then the general position condition is satisfied because for any pi G Pi and p2 G P2 
the line segment [pi,P2] is not parallel to a line segment contained in the unit sphere 
of this normed space (see Figure 3). In addition, d{Pi,P2) > 5 > 0. However, 
given e > 0, there can be no A > such that for any pair (P{,P2) of nonempty 
subsets, the inequalities D{Pi,P{) < A, D{P2,P2) < A imply D{Ri, R[) < e and 
D{R2, R'2) < e. Indeed, assume by way of contradiction that such A > exists. 
Let no G N be large enough such that 2/n < A for each n > tiq, n G N. Let 
P[ = {(0, lOn) : n< no, n G N} U {(2/n, lOn) : n > no,n G N} and P^ = P2. Then 
D{Pi,Pi) < 2/no < A and = D{P2,P^) < A. However, D{R[,Ri) = 00 since the 
strip S = {{xi,X2) : Xi < —10, X2 > lOno} is contained in Pi = dom(Pi,P2) but its 
intersection with R[ = dom(P{, P2) is empty. 

The example described above may seem somewhat complicated, but attempts to 
construct simpler ones may be futile: for instance, the simple setting where {X, d) — 
(M^ £00), Pi = {(0, 1)}, Pi = {(/3, 1}, P2 = {(0, -1)}, P^ = P2, and /3 > is arbitrarily 
small, imply that D{Ri,R[) = 00 (because the negative part of the horizontal axis 
belongs to Pi and not to P2) and there is no stabihty. However, in this case Pi and 
P2 are on a line segment parallel to a line segment contained in the unit sphere of the 
space. 

7. Proof of the stability theorem 

This section is devoted for proving Theorem 5.1. In the sequel, unless otherwise 
stated, (X, | ■ |) is a M™, m G N, with some given norm | • |; X is a nonempty compact 
and convex subset of X having more than one point; P,P',A, and A' are nonempty 
compact subsets of X; {Pk)keK and {Pl)k&K are two finite tuples of nonempty compact 
subsets of X representing the sites and the perturbed ones respectively; for each k E K, 
we let Ak = Ujyfc ^'k — Ujyfc Pjl unit vectors will usually be denoted by 9 or 

(f). The unit sphere of the normed space is denoted by Sj^. The distance between a 
set (or a point) and the empty set is defined to be infinity. In general, we note that 
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some of the claims and definitions given in the sequel hold in a more general setting 
(e.g., in some of them compactness or finite dimension arc not needed). 

The following 2 claims will be used along the way. They are special cases of claims 
formulated and proved in [49] and [50]. 

Theorem 7.1. (The representation theorem) Suppose that for all x E X the 

distance between x and P is attained. Then dom(P, A) is a union of line segments 
starting at the points of P. More precisely, given p & P and a unit vector 6, let 

T{e,p) = supO e [0,oo) -.p + te e X and d{p + te,p) < d{p + te,A)}. (2) 

Then 

dom{P,A)=[j [j\p,p + T{e,p)e]. 

peP 16*1=1 

Proposition 7.2. Suppose that mi{d(Pk, Pj) : j,k G K, j ^ k} > 0. Let e > be 
such that e < mf{d{Pk, Pj) : j,k E K, j ^ k}/Q. Suppose that the following conditions 
hold: 

3A e (0, e) such that for each k E K,p E Pk, y G dom{Pk, A^), andx G [p,y] 

if d{x,y) = e/2, then d{x,p) < d{x,Ak) — A. (3) 

3 A' e (0, e) such that for each k e K,p' e Pj,, y' e dom(Pj^, A'^), and x' e \p' , y'] 

if d{x', y') = e/2, then d{x' ,p') < d{x', A'^) - A'. (4) 

Suppose that there are M, M' e (0, oo) such that for all k G K, 

snp{Tk{e,p) -.pe P,\e\ = l} < M, snp{U{e',p') : p' G P', \e'\ = 1} < M', (5) 

where Tk{6/p) and Tl{6\p') are defined as in {2), but with A^ and A'f, instead of A. 
Let €4 be a positive number satisfying 

4(1 + M/e)e4 < A/2, 4(1 + M'/e)e4 < X' /2. (6) 

Let 

Rk = dom(Pfe, Ak), R'k = dom(P^, ^'J. 

D{Pk,Pl)<e4 ykeK, (7) 
then D{Rk, R'f^ < e for each k G K. 

Definition 7.3. Given two nonempty subsets P and A of X , the notation [P, A] is 
the set of all non- degenerate line segments of the form [p, a], where p E P , a E A. The 

notation [P,A] is the set of all unit vectors generated by segments from [P,A\. The 
notation P, A is the set of all unit vectors generated by endpoints of segments from 
[P, A] . In other words, 

\P^]-[^,--\p.a]E[P,A]\. (8) 
PrA = \P^]U[AJ^]. (9) 
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The notation A is the set of unit vectors generated from nondegenerate line segments 
contained in A, namely 

A= l-^L^ .\a,a']CA,a^a'\ . (10) 

The most important case for us for (10) is when A = Sj^. If for instance the 
normed space is (]R^,£oo)) then = {(0, 1), (0, — 1), (1, 0), (— 1, 0)}. In this case Sj^ 
is compact, but in general it may not be, and this is the case of from Example 3.3. 
If the unit sphere does not contain any line segment (i.e., the normed space is strictly 
convex), then = 

The following lemma establishes simple properties of P, ^. 

Lemma 7.4. (I) Given x,y & X , x ^ y, the segment [x,y\ is parallel to some 
segment contained in [P,A\ if and only if {y — x)/\y — x\ E P,A, i.e., if and 
onlyif±{y-x)/\y-x\e[P,A]. ^ 
(II) Let H = {t{a - p) : a e e P, i e R}. Then P,A^HnS^. 

Proof. Part (I) follows directly from the fact that the non-degenerate segments [x, y] 
and [p, a] are parallel if and only if {y — x)/\y — x\ — {a—p)/\a—p\ or —{y — x)/\y — x\ — 
(a—p)/\a — p\. Part (II) is simple too: if = ±{a — p)/\a — p\ G P,A, then |^^| = 1 and 
e = t{a-p) for t = ±l/\a-p\. Hence 9 e HnS^. On the other hand, if 9 e HnS^, 
then 9 = t{a — p) for a G A, p G P, and t G M. Since |^^| = 1 we have a ^ p and 
\t\ = l/\a — p\. Hence either t = l/\a — p\ or t = —l/\a — p\ and thus 9 & P,A. □ 

The following lemma establishes a simple estimate on the Hausdorff distance of 
unions of sets in term of the Hausdorff distances of members in the unions. 

Lemma 7.5. Let {Gi)i^j and {G'^i^i he two tuples of subsets of a metric space, where 
I is any nonempty index set. Then 

D{U,eiG,, U,eiG',) < sup{D{G,, G[) : t e I}. 

Proof. Let y G Lli^jGi. Then y G Gj for some j G /. Hence 

d{y, UieiG'^ < d{y, G'^) < sup{d{x, G'^) : x G Gj} < D{Gj, G'^) < sup{P>(G',, G'^ : i e 1} 

by the definition of the Hausdorff distance D. Thus sup{d{x, Uig/G-) : x G Uig/Gj} < 

sup{D{Gi,G'j) : i E I}. In the same way sup{(i(a;', Ujg/Gj) : x' G Ujg/G^} < 
sup{P'(G'i,G'^) : t e I}. Therefore D{Ui^jGi,UiaG'^ < sup{D{Gi,G'^ ■ i e I}, 
again by the definition of the Hausdorff distance. □ 

The next lemma establishes estimates on \d{P',A') — d{P,A)\ and D(P, A, P', A') 
whenever there arc known estimates on D{P,P') and D{A,A'). 

Lemma 7.6. Suppose that D{P, P') < ei and I?(A, A') < 62 for some positive numbers 
61,62. Then 

(I) \d{P\A')-d{P,A)\< 61 + 62 
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(II) Ifei + 62 < d{P,A) then 



and 



D{[P,A],[P',A'])< 



2(ei + 62) 
d{P, A) 

2(61 + 62) 



(11) 



(12) 



d{P, A) 

Proof. Wc first prove part (I). Let p' e P' and a' e A' be given. Since D{P, P') < 61 
and D{A, A') < 62 there are p & P and a & A such that d{p,p') < 61 and rf(a', a) < 62, 
so 

d{P, A) < dip, a) < d{p,p') + d{p', a') + d{a', a) < 61 + 62 + d{p', a'), 

Since p' and a' were arbitrary it follows that d{P, A) — 61 — 62 < d{P', A'). In the same 

way d{P', A') - 61 - 62 < d{P, A). 

We now pass to part (II). It suffices to prove (11) since in the same way an analogous 

inequality with [A.-P], [A',P'] is proved, and using Lemma 7.5 and P,A = [P,A\ U 

fVP] we obtain (12). 

Since 61 + 62 < d{P, A) it follows that \a—p\ > and, from Part (I), that \a' —p'\ > 
whenever p E P, a E A, p & P' , a' E A' and d{p,p') < 61, d{a, a') < 62. Now suppose 
that [p,a] G [P,A] is given. Since D{P,P') < 61 there exists p' G P' such that 
d{p,p') < ei- Since D{A,A') < 62 there exists a' G A' such that d{a,a') < 62. Hence 



P' 



p + hi(pi and a' — a + /i202 where hi G (0,6j) and = 1, 



1,2 (e.g.. 



01 = (p' — p)/\p' — p\). Let h — /i202 — ^101 and u — a — p. Then a' — p' — u + h and 



a' — p' 


a — p 




\a — p (a' — p') — \a' — p'\{a — p) 


|a' — p' 


\a — p\ 




\a' — p'\\a — p\ 



u\{ 


u + h)- 


- u 


u + h\ 


_ \{u + h){\u 




u + h\) + 


u + h\{u + h — u)\ 




a' — p'\ 


\a - 


-P 






\u + h\ 


\a - 


-P 





< 



2\u + h\\h\ 2(61 + 62) 
\u + h\d{P,A) d{P,A) ^ ' 

where the triangle inequality was used a few times for establishing inequality. In the 
same way as in (13) for any [p\ a'] G [P', A'] there exists [p, a] G [P, A\ such that 

(ei + 62) 



a — p a — p 



\a' — p' 



\a — p\ 



< 2 



d{P, A) 



Therefore we obtain (11). 



□ 



The next lemma proves several equivalent conditions for having equality in the 
triangle inequality. This lemma will play an important rule later. It is based on the 
following recent theorem of L. Maligranda. 

Theorem 7.7. (iMaligranda [44, Corollary 1.5]) Let Xi and X2 be two non-vanishing 



vectors in a nornied space. Let Xi 
and only if \xi+ — 2. 



1,2. Then \xi + X2\ = \xi\ + \x2\ if 
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Lemma 7.8. Let xi and X2 be two non-vanishing points in a normed space. Let 
Xi = Xi/\xi\, i = 1,2. Then the following conditions are equivalent: 

(a) The segment [xi,X2] is contained in the unit sphere. 

(h) \xi + X2\ = \xi\ + \X2\. 

(c) \xi +^2! = |^i| + 1^1- 

(d) \y\ — 1 for some nontrivial convex combination y of xi and £2. 

Proof, li Xi = £2, then all the conditions hold trivially. Assume from now on that 
Xi 7^ X2- If is contained in the unit sphere, then so is the convex combination 

of the form 0.5xi + 0.5x2- Therefore \xi + = 2, and hence \xi + X2\ = + \x2\ 
by Theorem 7.7. Hence (a)^»(b). The assertion (b)=4»(c) follows from Theorem 
7.7. The assertion (c)=^(d) is trivial because the equality |fi +^2! = |^i| + |^| is 
equivalent to |0.5a^i + 0.5a52| = 1. 

It remains to prove the assertion (d)^^(a). Let y be a nontrivial convex combina- 
tion of Xi and £2 satisfying \y\ = 1. Suppose by way of contradiction that some point 
z in the nonempty segment {xi,y) is not on the unit sphere. Since the unit ball is 
convex, the segment [xi,y] is contained in the unit ball. Hence \z\ < 1. Since X\^ £2 
by assumption and since £\ 7^ —£2 (otherwise \y\ < 1 as can be easily checked) the 
real linear combinations of £1 and £2 span a (real) two dimensional vector space. 

Since y e {z, £2) we can write y = l3z + {1 — (5)£2 for some /? e (0, 1). Consider the 
point of intersection w between the ray {ty : i > 0} and the segment a!2], where 
z = z/\z\. It can be written asw = ty for some t > 0. From the linear independence of 
z and £2 and by equating coefficients it can be easily verified that t = 1 / {(3\z \ + 1 — (3) . 
Thus t > 1 and \w\ = t\y\ = t > 1. But w is a convex combination of the points z and 
£2 which belong to the unit ball. Hence \w\ < 1, a contradiction. This contradiction 
proves that any point in the segment {£i,y) belongs to the unit sphere. In the same 
way any point in (^,^2) belongs to the unit sphere. Hence the whole segment [xi,f2] 
is contained in the unit sphere as claimed. □ 

Remark 7.9. By appropriate modifications, e.g., considering the convex hull gener- 
ated by the vectors and checking whether it is contained in the unit sphere, the result 
described in Lemma 7.8 can be generalized to equality involving any finite number 
of non- vanishing vectors, and even infinitely many (the latter case in a Banach space 
context). 

The following lemma generalizes is based on the previous one and is needed for later 
use. 

Lemma 7.10. Let p and y be points in a normed space and let A be a nonempty 
subset contained in this space. Suppose that p ^ A and that d{x,A) is attained for any 
X G \p,y). Suppose also that d{y,p) < d{y,A) and that for each a & A the segment 
\p,a] is not parallel to any non- degenerate segment contained in the unit sphere of the 
space. Then the inequality d{x,p) < d{x,A) holds for any x & \p,y)- 

Proof. If y = p, then the assertion is obvious (void). Now assume that y ^ p and let 
X G [p,y). U X = p, then obviously d{x,p) < d{x,A) because p ^ A and d{p,A) is 
attained. Assume that x G {p, y) and let a G A satisfy d{x, a) = d{x, A). Then y a 
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because otherwise \y — p\ < d{y,A) < \y — a\ — 0, i.e., y — p, a. contradiction. In 
addition, since x & \p,y], 

\y — x\ + \x — p\ — \y — p\ < \y — a\ < \y — x\ + \x — a\. (14) 

Therefore |x— p| < |x — a|. Assume by way of contradiction that 

\x — p\ — \x — a\. (15) 

Then (14) imphes that |y — a| = \y — p\ and 

W ~ y\ = |a — a;| + |a; — |/|. (16) 

It must be that a is not on the ray emanating from y and passing via p. Indeed, 
in order to be on this ray there are two possibilities. If a is to the left of p, then 
|y — a| > \y — p\, a contradiction to \y — a\ = \y — p\; if a G {p,y), then \a — y\ < 
\y — p\, a contradiction to d{y,p) < d{y,A). Finally, from previous discussion we 
know that a ^ p and a ^ y. Thus a is indeed not on the ray and in particular 
a ^ X. Hence (a — x)/\a — x\ {p — y)/\p — y\ = {p — x)/\p — x\. It follows that the 
segment [{p — x)/\p — x\, (a — x)/\a — x|] is non-degenerate and by (16), the equality 
{x — y)/\x — y\ = {p — x)/\p — and Lemma 7.8 this segment is contained in the unit 
sphere. But this segment is parallel to [p, a] since (15) implies that 

a — X p — X {a — p) 

\a — x\ \p — x\ |p — x| ' 

and so 

{{a — x) /\a — x\) — {{p — x) /\p — x\) {o, — p) /\p — x\ a — p 

\{{a — x)/\a — x\) — {{p — x)/\p — x\)\ |(a — p)/|j9 — x|| |a — p| 

We arrived to a contraction since it was assumed that [p, a] is not parallel to any non- 
degenerate segment contained in the unit sphere. Thus — p| < — a| = d{x, A). □ 

The following lemma proves an estimate based on compactness. 

Lemma 7.11. Suppose that d{P,A) > 0, that A is compact, and that the following 
condition holds: 

for any p & P, a & A and x,y & X, 

if d{y,p) < d{y,a) and x & [p, y), then d{x,p) < d{x,a), (17) 

Suppose that e > satisfies e < d{P,A)/2. Then the following condition holds: 

There exists A e (0, e) such that for each p E P,y E dom(P, A), x E \p,y], 

if d{x, y) = e/2, then d{x,p) < d{x, A) - A. (18) 

Proof. Define 

S — {{x,p) : X E X,p E P and there exists y E X such that 

xE\p,yld{y,p)<d{y,A) &ndd{x,y) = e/2], (19) 

A = inf {(i(a;. A) — d{x,p) : {x,p) E S}. 
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Given p & P, let b E A he arbitrary. Let 6 = {b — p)/\b — p\, x = p + 0.25e6, and 
y = p + 0.75e9. By the convexity of X we have x,y E [p, b] C X. Given a E A, since 
d{P,A) > 2e we have 2e < d{a,p) < d{a,y) + d{y,p) = d{a,y) + 0.75e. Because a was 
arbitrary we have 1.25e < d{y,A). It follows that d{y,p) < d{y,A) and {x,p) e S. 
Thus 5 0. 

For proving the assertion it suffices to show that A is attained at some (x, p) G S. 
Indeed, let y correspond to x and p in the definition of S. By assumption A is 
compact. Hence d{x, A) = d{x, a) for some a E A. Prom the property of y we have 
d{y,p) < d{y,A) < d{y,a). Thus, by (17) it follows that A = d{x,a) — d{x,p) > 0, 
and by taking 

A = min{e/2,A} (20) 

we see that (18) holds. 

To show that A is attained, consider the sequences (a;„)^-^ and {yn)'^=i such that 
{xn,Pn) G S and A = limn^^{d{xn, A) — d{xn,Pn)). For each n G N let |/„ correspond 
to {xn,Pn) in the definition of -S". We can write Xn — Pn + tn&n where \9n\ — 1 and 
tn E [0,oo). In fact, d{yn,Pn) > d{yn,Xn) = e/2 and also 9^ = {yn - Pn)/d{yn,Pn)- In 
addition, t„ < diam(X). Since P is compact, we can pass to a subsequence and obtain 
that t = lim^^oo^n; and p = lim/^ooPn; for some t E [0, diam(X)] and p E P. Since 
X is compact, then also y — limi^^ym. for some y E X and a subsequence (n^J^^ 
of positive integers, and then 9 — limj^oo^n;. for & — in ~p)/d{y,p). In addition, 
X — limi_^oo Xni. ioY X = p + t9. 

As a result, x E [p, y\. Finally, d{y, x) = e/2 and d{y,p) < d{y, A) by the continuity 
of the distance functions. Thus {x,p) E S and A = d{x,A) — d{x,p) > 0. 

□ 

The next lemma establishes a simple property of subsets having a finite face de- 
composition. 

Lemma 7.12. // S has a finite face decomposition S = ul'tlFj,, then any non- 
degenerate line segment [s, s'] contained in S must be contained in for some i E 
{!,...,£}. 

Proof. First note that [s, s'] nF^+i = since otherwise the fact that F^+in(uf^iFj) = 
implies (using the fact that Fi is closed for each i) that a small part of [s, s'] around 
any s" E [s,s'] fl F^+i is contained in F^+i, a contradiction. Now consider s: there 
exists a maximal nonempty subset / C {1, . . . ,£} such that s E Fj for any i E L 
If / = {1, . . . ,£}, then because s' E Fj for some j E {!,...,£} it follows that both 
s. s' E Fj and then, by convexity, [s, s'] C Fj. Otherwise / is strictly contained in 
{!,...,£}. 

It must be that some s" E {s, s'] belongs to Fj for some i E L Otherwise there 
exists a subsequence (s^)^j^ converging to s but satisfying s'^ ^ Fj for any i E I and 
n eN. Because {1, . . . ,£}\I is nonempty and finite it follows that infinitely many of 
these s'^ belong to Fj for some fixed j ^ /. But Fj is closed by assumption, so the 
limit s is in Fj too, a contradiction to the maximality of J. Therefore there exists 
some s" E Fi n {s,s'] for some i E I, and by convexity it follows that [s,s"] C Fj. 
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Hence [s, s'] C Fj by the maximality property of Fj (see Property (b) in Definition 



The foUowing lemma estabhshes a simple property of a set which is related to 
Definition 7.3. 

Lemma 7.13. Suppose that F is a nonempty convex subset of a real vector space V. 
Let 

H = {t{v — u) : t eM., u.v E F are arbitrary}. 

Then H is a linear subspace ofV. In particular, ifV is finite dimensional, then H is 
closed with respect to any norm defined on V . 

Proof. Let uq E F he given and let 

t 

Ui & F are arbitrary for each i} 



be the linear subspace spanned hy F — uo — {f — uq : f E F} . It suffices to prove that 
G = H. The inclusion H C G is simple, since t(v—u) = t{ui—UQ) + {—t){u2—uo) where 
Ml = V, U2 = u. Now consider the inclusion G ^ H. Given a finite linear combination 
w = Eti - uo) e G, let / = {1, . . . , i}, I. = {i e I : Xi < 0}, 1+ = {i e I : Xi > 
0}, /o = {i G / : Ai = 0} and a = 1/(1 + ^jgj |Aj|). For each i E I let Vi — Ui — Uq 
and Pi = a\Xi\. Finally, let u' := Xlig/_ A"^* '■~ Xlig/+ A"^* with the agreement 

that whenever the sum is over the empty set, the sum is considered as 0. With all 
of these notations it can be easily verified that u' + aw = v' (this is actually true 
independently of the specific value of a). As a result, w — {l/a){{v' + uq) — (it' + ito)), 
so for finishing the proof it suffices to show that u' and v' are in F — mq- 

Since F is convex, then so is F — uq. This implies that u' and v' belong to the 
set F — Uq. Indeed, consider for instance u' . Let 7 = « Ejg/_ |Ai|. If 7 = 0, then 
u' — Q E F — Uq. Otherwise < 7 < 1 and hence 7^1 E F — uq since this element 
is a convex combination of Vi and 0. Finally, since u' — X]jg/_('^l'^il/7)(7'^i) ^^<i 
1 = Ejg7_ (a|Aj|/7), it follows that u' is a convex combination of elements from 
F — Uq. Therefore u' E F — Uq. 

Finally, of in addition V is finite dimensional, then it is well-known that any linear 
subspace of V is closed with respect to any norm defined on V. □ 

The next few lemmas establish additional estimates needed in the proof of the main 
result or in the proof of auxiliary assertions. 

Lemma 7.14. Suppose that d{P, ^4) > and any segment [p, a] E [P, A] is not parallel 
to a line segment contained in the unit sphere of the space. Suppose also that 

has the finite face decomposition Sj^ = U^^jFj. Then d{P, A, S^) > 0. 

Proof. The assertion is obvious (void) if Sj^ = 0, since the distance to the empty 
set is infinity. Now assume that Sj^ % and suppose by way of contradiction that 
d{P., A., S^) = 0. Then there exist sequences (0„)5^i, {On)^=i such that lim^^oo |0„ — 



3.1). 



□ 




1=1 
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9n\ = 0, (pn ^ P,^, and 9n G for each n E N. By the definition of for each 
neN we have 6'„ = (s', - s„)/(s^ - Sn) where [s^, s'J ^ Sj^, Sn ^ s',. 

Since P,A= [P,A] U there exists an infinite subset Nq of N such that 

4>n £ [P, A] for each n G Nq, or 0„ G [A, P] for each n G Nq. Assume the first case. 
The other case can be considered in the same way. Then there exist sequence {pnjneNo 
and {an)neNo such that a„] G [P, A] for each n & Nq and 



hm 



I C^n I 



0. (21) 



By passing to subsequences and using the compactness of P and A we conclude that 
there exist p & P and a & A such that p — limn^oo,neNi Pn and a — \i'mn^oo,neNi o,n for 
an infinite subset A^i of A^o- By Lemma 7.12 for each n e Ni there exists i G {1, . . . , ^} 
such that [sn, s^] C Pj. Because there arc finitely many subsets Pj, i G {!,...,£}, there 
exist i G {!,...,£} and an infinite subset N2 of A^i such [s„, C Pj for each n G iV2. 

Consider the set H = {t{v — u) : t E M., u,v E Pj}. By Lemma 7.13 and since 
the space A" is a real finite dimensional the set H is closed. By Lemma 7.4(11) the 

set Pj, Pj is nothing but the intersection of H and the unit sphere S^. Hence Pj, Pj 
is closed. Using (21) with n E N2 we conclude that {a — p)/\a — p\ is a hmit point of 
elements from Pj, Pj. Thus (a — p)/\a — p\ E Fi, Fi. Hence {a — p)/\(i — p\ = t{u — v) 
where u,v E Fi,u 7^ v and t G M. Either t = —l/\u — f | or t = l/\u — v\ and hence 
(a — p)/\a — p\ = ±{v — u)/\v — u\. Because P^ is convex the line segment [u,v] is 
contained in Pj C 5"^. Prom the preceding sentences it follows that [p, a] is parallel 
to [u, v] which is a nondegenerate line segment contained in the unit sphere. This is 
a contradiction to the initial assumption. Hence d{P,A,Sj^)>0. □ 

Lemma 7.15. Suppose that mm{d{Pk, Pj) : k,j E K,k j} > and that no two 

points from different sites form a segment which is parallel to a non- degenerate segment 
contained in the unit sphere of the normed space X . Then there exists r > such that 
d{Pk,Ak, Sj^) > r for all k E K where A^ = Uj^^Pj- 

Proof. By Lemma 7.14 we know that r^j := d{Pk,Pj,S^) > for any k,j E K, 
j 7^ k. (If Sj^ — 0, then any positive number, say 1, satisfies the assertion.) Hence 
minjrfcj : k,jE K, j ^ k} > 0. Now Hx k E K and let (j) E Pk, A^ and 6 E S^^. Then 
(f) = {pj - Pk)/\Pj - Pk\ for some pk E Pt and pj G P,-, or = (g^ - qj)/\qk- Qjl for 
some qk E Pk and qj E Pj. Thus |0 — > t^j. Hence d{Pk. Ak, Sj^) > minlr^j : j 7^ 
k,j E K}. Therefore any positive r satisfying r < min{rA;j : j,k E K,j ^ A;} satisfies 
the desired conclusion. □ 

Lemma 7.16. Let e > and r > be given. Let 

^€,r — {{P:"^) '■ P a-iid ^are nonempty compact subsets oiX satisfying 

d{P, A)/2>e and d{P^, %) > r}. (22) 
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Then there exists a universal constant X e (0, e) such that (18) holds with this A for 
any {P,A) G ^e,r- 

Proof. Note that the A described in Lemma 7.11 depends on P and A, while here 
it is independent of them. To avoid ambiguity, we denote the A from Lemma 7.11 
by Xp^A- Given any {P,A) e ^e,r it follows by Lemma 7.10 that (17) holds since 
d{P,A, S^) > r. Hence Lemma 7.11 does imply the existence of Xp^a > such that 
(18) holds. Recall also that Xp^a < e- Define 

X^mi{Xp,A:iP,A)e^,^r}- (23) 

For finishing the proof it suffices to prove that A > 0. Suppose by way of contradiction 
that A = 0. Then for any n E N there exist compact subsets P*^"-* and A^^^ of 

X satisfying <i(P("), A("))/2 > e and d(P(^^rA("), %) > r such that the number 
A„ := Ap(n) A(") > corresponding to them from Lemma 7.11 satisfies A„ < 1/n. Note 

that by (19) and (20) we have A^j = A„ < e/2 for large enough n. This, (19), and 
the compactness of A^") imply that there exist p„ e P("\ y„ E X, Xn E [pn,2/n], and 
a„ e such that d{xn,yn) = e/2, d{yn,Pn) < (^(z/n,^^"^) = d{yn,an) and 

An < d{xn,an) -d{xn,Pn) < A„ + l/n < 2/n. (24) 
Because d(P("), %) > r we also have 

\{{an- Pn)/\an- Pn\) " ((^2 " "5l)/|s2 " Sl|)| > T (25) 

for any nondegenerate segment [si,S2] C Sj^. By passing to convergent subsequences 
we obtain the existence of points p,x,y, and a contained in X such that x E [p,y], 
d{x,y) — e/2, d{y,p) < d{y,a), d{p,a)/2 > e and d{x,a) — d{x,p) = 0. If x = p, 
then p = x = a E A, contradicting the fact that d{P, A) > 0. Hence d{x, a) = 
d{x,p) > 0. Exactly as in Lemma 7.10 (from (14) until the end) it follows that [p, a] 
is parallel to a nondegenerate segment [si,S2] that is contained in the unit sphere 
Sj^ (the hue segment [{p — x)/\p — x\, (a — x)/\a — x\]). But this is impossible since 
|((a —p)/\a — p\) — {{s2 — •Si)/|s2 — si|)| > r because of (25), a contradiction. □ 

We are now ready to prove Theorem 5.1. 

Proof of Theorem 5.1: It is sufficient to assume that e G {0,r]/6), otherwise one 
can take for the given e the A associated with, say, 77/ 10. We will show that all the 
conditions needed in Proposition 7.2 are satisfied. 

Since X is bounded we obtain (5) with, say, M — M' — diamX. Let k E K he 
given. Since no two points of Pfc and Ak — VJj^kPj form a segment which is parallel 
to a non-degenerate segment contained in the unit sphere, Lemma 7.15 implies the 
existence of r > such that d{Pki A^, S^sj^) > 2r for each k E K. By (1) it follows 
that d{Pk,Ak) > 77 > 2e for each k E K. Therefore {Pk,Ak) E '^e,r (see (22); note 
that we take here ^^^^ and not ^'e,2r) for each k E K. 

Let A be taken from Lemma 7.16 and let A' = A. Let A be any positive number 
satisfying 

A • f ^ ^^1 

A < mm < — ; — > . 

\8(l + (M/e))' 8 / 
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In particular A < A < e. Let {Pl?)k<^K be any tuple of nonempty compact subsets 
of X satisfying D{Pk,Pl^) < A for any k G K. For each k E K let A'j. = Uj^kPj- 
Then D{Ak,A'^) < mm{D{Pj,P^) : j e K, j k} < A < e hy Lemma 7.5. This and 
Lemma 7.6(1) imply that 

diPl, A[) > diPk, Ak)-2e>7]-2e> 2e. (26) 

In addition, Lemma 7.6(11) implies that D(Pk, Aj., Pl^, A'jJ < AA/rj < r. 

Fix k G K and let s G S^m and q' G P^, A'f^ be given. From the inequality just 
proved, namely D{Pk, Ak, PI., A'jJ < r, there exists q G Pk, A^ such that d{q,q') < r. 
In addition, d{q,s) > 2r because d{Pk, Ak, Sjs^m) > 2r. Hence 

2r < d{q, s) < d{q, q') + d{q', s) < r + d{q', s). 

As a result, d{PI, A'^, Sum) > r since q' and s were arbitrary. Consequently (see (22)), 
this and (26) imply that (P^, A'^.) G for each k G K. Hence Lemma 7.16 implies 
that (18) holds with P = P^ and A = A^, and also with P = Pj^ and A = A'j^ for each 
k e K. Therefore (3) and (4) hold with A' = A. Thus, by denoting 64 = A we see that 
all the assumptions in Proposition 7.2 are satisfied, and hence D{Rf;, R'^)) < e holds 
for each k e K. □ 

8. Concluding Remarks 

The stability property established in Theorem 5.1 is quite general in the sense that 
the norm is not restricted to be uniformly convex and the sites can have a pretty 
general form (they only need to be compact). However, in some aspects this result 
is quite restricted comparing to the one presented in [50] since there infinitely many 
sites were allowed, the sites and the world X were not assumed to be compact or 
bounded (but a certain boundedness condition was assumed), and the normed space 
could be even infinite dimensional. Moreover, explicit (dimension free) estimates 
were given there in contrast with the non explicit ones given here. In particular, 
it is not clear whether A given in Theorem 5.1 does not depend on the dimension. 
However, we feel that by strengthening our approach it is possible to achieve the above 
properties in more general normed spaces, but for this one has to find ways to obtain 
explicit estimates not based on compactness or finiteness arguments. In addition, we 
also feel that the "general position" assumption on the sites can be weakened, with 
some caution. It will be interesting to establish such results. In particular, it will 
be interesting to weaken the finite face decomposition assumption imposed on the 
structure of the unit sphere and to see whether a stability result can be formulated 
for the normed spaces induced by the unit spheres mentioned in Examples 3.4-3.5. 

It is also interesting to generalize the results to other settings, such as manifolds, 
weighted distances, and convex distance functions. Finding theoretical and real-world 
applications of the result (in addition to the possible ones described in Section 1) may 
be of interest too; a promising place where such a result can be applied is in the con- 
text of non-Euclidean stochastic geometry and other kind of probabilistic questions 
related to geometry and the distribution of the sites, as done in the Euclidean case 
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(see, e.g., [45, pp. 39, 291-410] for a discussion related to Poisson process). Finally, it 
may also be interesting to discuss the possibility of stability where there is no one-to- 
one correspondence between the original sites and the perturbed ones, e.g., because 
there were merges or eliminations due to some processes. A corresponding situation 
occurs when considering point clouds and it may have applications in data analysis 
and reconstruction, as can be seen in the somewhat related discussion in [15] (in the 
context of Euclidean lambda-medial axis). 
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